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Abstract. We consider the asymmetric simple exclusion process in d > 3 with open 
boundaries. The particle reservoirs of constant densities are modeled by birth and death 
processes at the boundary. We prove that, if the initial density and the densities of the 
boundary reservoirs differ for order of e from 1/2, the density empirical field, rescaled 
| as converges to the solution of the initial-boundary value problem for the viscous 

QJ ■ Burgers equation in a finite domain with given density on the boundary. 

oo ; 
> ■ 

t^j- ■ 1. Introduction 

A driven lattice gas with open boundaries is a system of particles jumping at random on 
a lattice, subject to the action of an external field and exchanging matter with a reservoir 
■ at his boundary. The combined action of the force field and the density gradient induced 
by the boundary conditions forces the system to reach a stationary non-equilibrium state. 
Systems of this kind show a complex behavior exhibiting non-equilibrium phase transitions 
Ph| [K], [SZ]. The simplest example of driven lattice gas is the asymmetric simple exclusion 

43 ' (ASEP). In the one-dimensional totally asymmetric case an explicit stationary solution 

is known [D] showing a very reach phase diagram with different behavior of the steady 
current depending on the values of the fixed densities on the boundaries. 

In this paper we study the time-dependent measure of the ASEP with open boundaries 
in d > 3 in the macroscopic limit. The system is contained in a finite cylinder A e = 
[— e -1 ^ -1 ] x 7Tg _1 , with vr^ _1 the (d — l)-dimensional microscopic torus of size 2e~ 1 + 1 
with the axes in the direction x±, namely we impose periodic boundary conditions in all the 
directions but X\. In the bulk particles jump to one of the nearest neighbors if empty with 
jump rate p ei (p_ e J in the direction (— ej). We assume that the vector of components 
St = p £i — p_ ei is such that 5i > 0. On the boundaries X\ = —e^ 1 ^^ 1 we allow for 
production and destruction of particles in the following way. Let b(u) be smooth functions 
on [—1, 1] x nf^ 1 . A particle is added independently in each site of x\ = — e -1 , when the site 
is empty, with rate <5i(l/2 + eb(ex)) and removed independently in each site of x\ = e^ 1 , 
when the site is occupied, with rate 5i(l/2 — eb(ex)). If b{u) has constant values b±, 
this choice of the jump rates corresponds to coupling the system to reservoirs of constant 
densities l/2 + e6_ and l/2 + eb + respectively. In the paper we write the computations for 
this cylinder geometry, but the proof could be extended to a more general convex domain. 
The initial measure is the local equilibrium corresponding to a density profile which is a 
perturbation of order e of a constant profile: po(ex) = 1/2 + em (ex). 
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We prove a law of large numbers stating the weak convergence of the rescaled empirical 
field (r)t(x) = 0, 1 is the occupation number in the site x at time t) 



e^ 1 -±)*(es) 

xeA £ 



with S(-) the Dirac measure, to the solution of the viscous Burgers equation in the domain 
Ai with density b(u) on the boundary Y 

d t m(t, u) =5- W u (m(t, u)) 2 + Y?i, j= i D i,jdl t)U m{t, u) 
777,(0, •) = mo(-) 
m(t, • ) = b{ • ) for t > 0, 

where D is a positive definite diffusion matrix whose expression is given by the Green- 
Kubo formula for this model. We remind here that the transport matrix is supposed to be 
infinite in d < 2 [Sp]. 

This special choice of the initial condition as well as of the boundary conditions (the 
difference between the top and the bottom densities is of order e) is forced by the fact that 
we want to study the behavior of the system on the diffusive time scale e~ 2 to see the effect 
of a finite dissipation. On this time scale the transport term, which involves first order 
time derivatives, is enhanced by a factor of order e~ l so that if it is of order e at time 
its contribution stays finite in the limit. This is also called incompressible limit [EMY1] in 
analogy with the Navier-Stokes case [EMY2] . 

There are few rigorous results on the hydrodynamic limit for interacting particle systems 
in a bounded domain. In [ELS] it is proved the hydrodynamic limit for one-dimensional 
gradient systems both in the time-dependent and stationary case. In [KLO] the analogous 
result has been obtained in the stationary case for one-dimensional non-gradient models. 
Finally, in [LMS] the latter result has been extended to the (i-dimensional case. All these 
papers deal with bulk reversible dynamics. If the lattice gas is driven by the boundary 
conditions, in general the stationary measure does not coincide with the invariant measure 
of the dynamics in the infinite volume, so that even if the generator of the bulk dynamics is 
reversible versus its invariant measure, the total dynamics is not. The asymmetric simple 
exclusion is an example of a lattice gas driven by an external field and it is not reversible 
w.r.t. his invariant measure in the infinite volume, which is a product measure. The 
case we consider in this paper is an example of a lattice gas driven by both external field 
and boundary conditions, so that there are two sources of non-reversibility. The effect 
of the asymmetry is seen in the macroscopic current as the term 5im 2 which gives rise 
to the transport term in the macroscopic equation. We have to face the difficulty of non 
reversibility in the bulk and the fact that the stationary measure is not explicitly known. 

The first problem has been solved in [EMY1] where it is proved the law of large numbers 
in the incompressible limit for the system without reservoirs. The method used there is 
based on the use of relative entropy of the true measure with respect to some suitable local 
equilibrium measure and the main point is to prove that this relative entropy vanishes in 
the limit e — > 0. In this case it is natural to assume for the local equilibrium measure 
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a product measure because this is the invariant measure of the generator. However, for 
the open system it is not clear which is the good candidate for describing the system for 
e small. We introduce a measure fi £ which is a product measure with chemical potential 
e\(t,ex) such that E^ £ [r] x ] = 1/2 + em(t, x), with m solution of the macroscopic equation. 
We prove that the relative entropy s{ji\ji £ ) of the non-equilibrium measure /x w.r.t. ji £ 
satisfy 



which implies the law of large numbers. One can understand the factor s~ 2 noting that the 
specific entropy of /i £ is of order e 2 . To get this result, we introduce an auxiliary measure 
jl £ which differs from \i £ near the boundary and prove the previous limit for the relative 
entropy w.r.t fi £ . Then, the result easily follows from the fact that lim^o e~ 2 s(/i £ \jl £ ) = 0. 

The method used in proving the hydrodynamic limit in the papers quoted above is 
different from ours because it is based on the martingale approach. However, in both 
methods it is needed an a priori estimate of the entropy (which assures also the control 
of the Dirichlet form) of the state of the process with respect to the invariant state of the 
full dynamics. Since this is not known, one uses some trial reference measure which is not 
invariant and as a result this entropy does not decrease in time. The difficulty to get an 
entropy bound in the asymmetric case is due to a diverging contribution to the flux of 
entropy due to the asymmetric part of the generator. This flux has to be controlled by 
the boundary generator. The choice of this generator is different from the one used in the 
case of reversible bulk dynamics, where the boundary death and birth process is required 
to satisfy a detailed balance condition. Our generator is a generalization of the one used 
by Derrida [D]. Our method would also provide the results for a reversible boundary 
generator, but then we have to speed it up much more than e~ 2 . 

One technical remark: the control of the terms in the bulk is done by using the replace- 
ment Lemma and the non gradient result in [EMY1] . This requires a localization procedure 
for the currents, which generates, in an open system, boundary terms. These terms are 
dealt with, in a reversible non-gradient case [LMS], by a simple use of integration by part 
and Schwartz inequality. The version of the integration by parts for the asymmetric case 
(Lemma 6.1 in [EMY1]) involves a function which is not bounded but just summable and 
as a consequence we need rather to use a modified version of the replacement Lemma to 
control these boundary terms. 

We conclude this section by noticing that the extension to the stationary problem of 
the convergence result in the macroscopic limit could give some insight on the structure of 
the non-stationary states. Recent papers [DLS], [BDGJL], focus on the problem of char- 
acterizing the stationary non-equilibrium measures in terms of large deviation functional. 
Unfortunately, the relative entropy method does not seem to be useful for the stationary 
problem, because it relies on estimating the time derivative of the relative entropy in terms 
of the entropy production which is expressed back again in terms of the relative entropy, 
so that one gets a closed differential inequality for the relative entropy. We think that the 
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usual martingale approach could be more suited for the stationary case. We plan to refer 
on that in a future paper. 

The plan of the paper is the following: in Section 2 we state the results. Section 3 is 
devoted to the strategy of the proof of the entropy bounds and to the control of the bulk 
terms while in Section 4 we have collected the estimates of the boundary terms. Finally in 
Section 5 we discuss extensions to different geometries and different boundary conditions. 

2. Notations and results 

Let e > such that e" 1 is integer and A e = {— e^ 1 , . . . , x n^ 1 be the cylinder in 
Z d , d > 3, of length 2e^ + 1 with basis rcf' 1 , the (d — l)-dimensional microscopic torus of 
size 2e~ 1 + I. We denote by r e = {x e A e | x\ = ±£ -1 } the boundary of A e . The elements 
of A £ will be denoted by letters x,y, . . . 

A particle configuration is described as an element rj e X £ = {0, l} Ae , where r)(x) = 1,0 
means that the site x is occupied by a particle or is empty. In this paper we are interested 
in the asymmetric simple exclusion process (ASEP) on A e with stochastic reservoirs at the 
boundary of A e . This Markov process is defined through its infinitesimal generator 

C £ = C £fi + C e>b . (2.2) 

The generator £ £j0 is the nearest neighbor ASEP. Its action on functions / : X £ — > R is 

d 

£e, f(v) = Yl E r x , x+e M[f(v x ' x+ei )-m], 

1=1 x,x+ei£A s 

where (ei, ... , ej) is the canonical basis of lR d and the rate functions r xx+e .(rj) are given 
by 

r x ,x+ ei (v) =Pe l v( x )( l -V(x + ei)) + p- ei rj(x + ei)(l - rj(x)) 
and i] x ' y is the configuration obtained by exchanging the occupations of sites x and y: 

{r)(y) \tz = x, 
r)(x) iiz = y, 
r}(z) Mz^x,y. 

For 1 < i < d, p ei ip-a) is the jump probability of a particle in the direction (—ei). 
For convenience we normalize the p e 's so that, for % — 1, . . . , d, p ei + p_ ei = 2. We call 
$i '■= Pei — P-en i — 1, • • ■ , d, and we suppose that 5i is non vanishing. We assume it 
positive without loss of generality. 

The currents W Xji are defined as 

w x M = -wS + 

rj(x) + r](x + ei) 



:= - [i](x + - i](x)] + Si 



r](x + ei)r)(x) 



Let A be the set [—1, 1] x n d 1 , where n d 1 is the (d — l)-dimensional torus with length 
2, whose elements are denoted by u, v, . . . We will consider as reference laws the Bernoulli 
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product measures v p on X £ which are defined, for any smooth function < p < 1 on A, 
by i/ p {j](x) = 1) = p{ex). It is well known that, in infinite volume, the product measure 
of Bernoulli laws with any constant parameter is invariant for the ASEP, but this is no 
longer valid in finite volume. 

The generator C £ ^ in (2.2) is the infinitesimal generator of a birth and death process 
which creates particles on the "left" of A e and destroys them on the "right". For any 
function / on X £ 

(v) = J2 cb ( £x ^)[f(v x )-m], 

xer £ 

where r] x is the configuration obtained from 77 by flipping the occupation number at site x 



r) x (z) 



r)(z) if z 7^ x 
1 — t](x) if z = x. 



The rate functions are designed to fix the value of the particle density at the boundary of 
A e . Let T £i+ (resp. r e _) be the "left" (resp. "right") side of T e : 

r e ,+ = { x 2 ,...,x d ) , (x 2 , . . . , x d ) e tt^ 1 } , 

r e _ = {(e~ 1 ,X 2 , ...,X d ) , (%2, ...,x d )e TT^ 1 } . 

For a smooth function b on the boundary T of A that we fix from now on, we choose 
C b (ex,i]) such that a particle is added on x G r gj+ , when the site is empty, with rate 
Si(l/2 + eb(ex)) and removed in r e _, when the site is occupied, with rate 5i(l/2 — eb(ex)): 

C\ex, r]) = Si (1/2 + eb(ex)) (1 - r](x)) l{x e T £:+ } + ^ (1/2 - eb{ex))r](x) l{x E . 

To understand this choice, remark that if we take 6 = then the density 1/2 is stationary 
for the full dynamics C e , therefore the leading coefficient 1/2 in the rate is necessary for 
the incompressible limit to make sense. The term eb will fix the value at the boundary of 
the density perturbation with respect to the stationary value 1/2. 

We will study the dynamics defined by the generator C e under the diffusive space-time 
scaling (ex, e~ 2 t). We denote by rj t the particle configuration at time t of the associated 
Markov process. When the process starts from the product measure p £ = Vi/2+em(.), m a 
smooth function on A such that m| r = b, its law is denoted by ]P^"\ 

As in [Y] and [EMY1], the proof is based on the study of the relative entropy. Given 
two measures H\ and p 2 on X e (with pi absolutely continuous w.r.t. p 2 ), the entropy of 
Pi with respect to p 2 is defined as 

H{p 1 \p 2 )= [ ^p-\og^p-dp 2 . 
J dp 2 dp 2 

We also introduce the relative entropy s(pi\p 2 ) := e d H(pi\p 2 ). 
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For simplicity we call v := vy 2 the reference measure. Notice that the measure \i e has 
a density \& with respect to v given by 

V(v) = \ ex P \ £¥( m (zx))v( x ) > "(v), 

where 

and Z is a normalization constant. Then the relative entropy of the measure fis (where / 
is some probability density) w.r.t. /i £ can be written as: 



The probability fi £ (t) is the law of i] t when the initial distribution of the process is /i e . Its 
density with respect to v is denoted by f t . We now consider the solution m(t,u) of the 
following partial differential equation 

d t m(t, u) = J2t =1 Sid Ui (m(t, u)) 2 + Ei<ij< d A,^. u .m(£, u) 
m(0, .) =m (2.4) 

m(t, .) = b(.) for i > 0, 
r 

where the diffusion coefficients D it j will be defined later on. 

We introduce the measure ^i/2+ £ m(t,.) whose density w.r.t. v is given by 



^t(v) = i ex P I Y] eip(m(t,ex))r)(x) \ , 



with Z t a normalization constant. 
Our main result is: 

Theorem 2.1. Let s(t) = s(f t \^ t ). Then, for any t > 0, 

\ime~ 2 s(t) = 0. 

e-+0 

Using standard arguments based on the entropy inequality (see e.g. [Y] and [EMY1]), we 
get as a corollary of the previous theorem the incompressible limit of ASEP with stochastic 
reservoirs: 

Theorem 2.2. Let ix e t be the empirical measure 

Then 7r| converges weakly in F b £ ' m -probability to the solution m(t,u) to (2.4)- 
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3. Bounds on the entropy 

The strategy in proving Theorem 2.1. is based on the study of the time evolution of the 
relative entropy. Here the difficulty relies on the fact that the invariant measure for the 
full dynamics is unknown. 

Actually, we will not work directly with the entropy s(t) but with some approximation 
of s(t) obtained by a modification of the density ^t- Let m be the solution of the equation 
(2.4) and 7 be a smooth function defined on a neighborhood of A such that <9 Ul 7 has a 
compact support included in ] — 1, l[x7r( d_1 ) with 7(-)| r = We introduce a new parameter 
< 9 < 1 which will go to after e, we denote by CIq the set 

n g -.= [-1 + e, 1 - e\ x 7r (d_1) (3.5) 



= 1 and 

^28 



and we choose a smooth function x 6 i n a neighborhood of A such that 
X e = outside IV Then we define 

p(t,u) = m(t,u) X e (u) + 7 («)(1 " /(«))■ 

Notice that p(t, .) is identically equal to m(t, .) in f^e and equal to 7 on Q \ fig. 
Recalling that ip was defined in (2.3), we denote by X(t,u) the function 

X(t,u) = <p(p(t,u)), 

then it is easy to check that, for u e ^20, 

d 

d t X(t, u) — ^2 hdui (A(i, u)) 2 - ^ A,< Uj A(t^)=o(5), (3.6) 

i=l l<i,j<d 

where e _1 o(£) is a bounded function. Moreover, with a suitable choice of x^j we ma Y 
assume that the derivative d Ui X , % — 1, . . . , d are bounded, uniformly in 6. 

Following [EMY1], we modify suitably the density \& t in the definition of the relative 
entropy s(t): for an integer £, we set k = £e~ 2 ^ d and A^ = {— k, • • • , k} d . The normalized 
indicator function uj is defined as 



UO[X) 



\k k \- l t{xek k }. 



Let {F n } be a family of local functions (i.e. which depend on the particle configuration 
only through a finite number of sites) on X £ and 

d 

®n(v) = Yl Yl {d Ui X)(t,ex)(T x F n *uj), 
i=i xeA 2 F k 



for x G Z d , where r. is the shift operator on X = {0, l} zd , A 2k = {—e 1 +2k,e l —2k}x7i" 



r d-i 

and * is the convolution product. We will work with n fixed and will take the limit n — > 00 
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after the limit e ^ 0. Therefore the index n will be omitted for sake of shortness. The 
modified density \&f is defined as in [EMY1] by 

*f (v) = i ex p I £ ( A * u ) (*> ^0*0 + e2 $fa) > > 

Z < UeA e J 

where Z^ 7 is the normalization constant. 
Moreover, we define 



and \& t the density 

^t(v) = — ex P S E e ( A * ^) (*' £X )v( x ) f • 
Z t Uga £ J 

We finally introduce the functional D , which are closely related to the Dirichlet 
forms. For any non negative function h and any measure \i on X £ , 

d 



i=l x,x+ei£A 



(3.7) 



Then we have 



Proposition 3.1. There exists some constant K > s«c/i that for every t > 0, 

lim Urn fim" lim" <jV 2 Sl (t) + ATe d ~ 4 / , *f + A>( V 7 ^, 1 = 0. 

6* ^0 n->oo >oo e^O I JO ^ ' J 

Remark. Since the Dirichlet forms D and Db are positive, Theorem 2.1 will follow from 
Proposition 3.1 if we can show that 

fim Em fimV 2 \s(t) - sAt)] < 0. (3.8) 

0->-(H->oo £->0 

Consider the density with respect to v given by 

' (3-9) 



Vt(v) = ex P j £ £A ^' ea 0»7(a0 [ 
* UgA £ J 



with the normalization constant. It is proved in Lemma 3.2 of [EMY1] that for any 
< 6 < 1 

Vme- 2 [8(f t \*l)- 8l (t)]=0. (3.10) 
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Moreover a straightforward generalization of the proof of this lemma together with the 
entropy inequality leads to the following estimate: for any a > 0, there exists a constant 
c(a) such that 

e^[ s (t)-s(f t \^l)] <ae- 2 s{f t \^])+c(a)9 (3.11) 

where we have used the fact that \(t,u) = <p(m(t,u)) for any t > and u e f^e- Then 
(3.10) and (3.11) imply (3.8). 

Proof of Proposition 3.1. A simple computation shows that 

£ ~ 2 f t s ^)- £d ~ 4 J ftC e hghfdu-e d - 2 J ^fftdu 

Since the function $(?y) depends on the configuration rj only through the variables 
{77(^1, • • • ,Xd) , \xi\ < (e^ 1 — k) + Sf}, with sp the support of F, then 

e d - 4 J f t C £ \ogh F du = e d ~ A J f t C £fi \oghfdu + e d - 4 J f t C £ , h \ogh t du. 

We use the basic inequality 

a log 7 < — (y/a — v 6) +a — 6 


for a and 6 positive numbers. It is applied with a — hf (rj x ' x+ei ) and b — hf (rj) to manage 
with £ £j0 and with a = h t (r] x ) and 6 = /^(r?) to manage with £ £i0 , We get 

d 



(3.12) 



£_2 ^ s i(*) < E(e,t)-e d - A D (y/hf,^u) 

+ e d - 4 J (C £tb h t )%du-^e d - 4 D b (^/h t , % v ) 
with £■(£, t) := i?i + E 2 , where 

/C* r f* fy F 

J 4f 

and £* denotes the adjoint operator of C £ $ in L 2 (z/) . After a change of variables the term 
Ei can be rewritten as 

d ^ r ~ (3.13) 

+ £d ~ 4 E E * / W^+ei)-^))(/*-*f)(^ 

■^El + El 
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where 

Using again a change of variables at the boundary, the second term of E\ is equal to 

Ei = e d -*Y<^ j Hx)-{r,{x))^) ni {ex)f t {r,)du 

xeTe (3.14) 
= e d ~ 4 ^25 1 / (r](x) - (r](x))~ t )n 1 {ex)f t {r])du , 

where n = (ni, 0, • • • , 0) is the outward unit normal vector to the boundary surface T, and 
by standard manipulations can be bounded by the Dirichlet form and the generator at the 
boundary: 

Lemma 3.1. For all a > 0, 

- d ~ 4 5> / W*) - <^)>*>(^)/ t (^ < 

™ r- "P J 

< -e 



e 

ser £ 

_d-4 



j {C £ , b h t )%dv + e d -^D b (^h t , %u) + e 3 k 4 ^ , 
for some positive constant C\. 

Proof. For e > 0, denote by C k eb the following generator at the boundary 
(£J l6 /)fa) = *i E C b + (ex, V ) [f( V x ) - m] 

+ *i £ C^ex^^f^-f^)], 



where 

&+{ex,ri) = (l/2 + e(\*u)(t,ex)) (1 - rj(x)) 
C b _(ex,r]) = (l/2-e(\*u)(t,ex)) T)(x) 

Remark that the replacement of C £tb by C k b in the derivative of the entropy produces a 
term that we can bound by the Dirichlet form. Indeed, we have 

e d ~ 4 ! ((C £yb -£ k £b )h t )y t dv 




= e d ~ 3 £ (b(ex)-(\*u)(t,ej) [(l-rtx^htirj^-hti^tdis 
+ e d ~ 3 J2 ((A *w)(M) J v(x)[h t ( V x ) - h t (rj)]%du , 
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by Taylor expansion at the second order, this last quantity is equal to 



e d - l k 2 G ^ ex ) [(l-r)(x))Mr] x )-ht(ri)}%dis 



</-4 a r^ i.rr ,7, ,.\ , ^3,4^1 



<£"-A(AM+« 

Zi CI 

for all a > and some constant Ci, where G(t,ex) is bounded and proportional to the 
second derivative of X(t, .). We have used, in the last inequality, Schwartz inequality and 
the fact that h t is a probability density with respect to \l/ t z/. 

To conclude the proof of the lemma, it is easy to check that the left hand side of the 

inequality in Lemma 3.1 can be rewritten as — e d ~ A f (c k e b h^ ty t dv , by using the change 
of variable formulas valid for any measure u g 

J {l-V*)f(v x ' + )dvM = 1 ~^ X) j V(x)mdu g (rj) 

j Vx f(r)*>-)dvM = j {1-V(x))f v du g ( v ) 

where i] x,+ and are the configurations obtained from 77 by adding or removing a particle 
at the site x. □ 

The term E\ is shown, by using Taylor expansion and Schwarz inequality, to be less or 
equal to 

/ (EE {F' i \*u)(t,ex)W Xti - £*$\ {f t -^()( V )dv + R £ (t) +o £ (l). (3.15) 

J I i=l x&A e J 

where W x>i is the current defined in section 2, c^A is the discrete gradient of A 

<9 e,i A(t, u) = e' 1 X(t, u + eet) - X(t, u) 

and for any function g on X e , we denote by g = g— < g > v . The operator Cq is the 
generator of the ASEP in infinite volume space X = Z x ti^ 1 and Cq is its adjoint in 
L 2 (z/), given for any local function / by 

In equation (3.15) as in the sequel, R £ (t) is real sequence of the form 

Rett) =e d J2 G ^ ex ) I r x g(v)(ft ~ *?)iri)dv (3.16) 

xeKe J 



£ d-2 
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with some bounded functions G(t, .) denned on A and g on X £ . Observe that, for such 
sequences, one can use the entropy inequality and obtain, for any A > 0, 

Rett) < Ae d - 2 Sl (t) + A-'o^l) , (3.17) 

for a real sequence o e (l) which is bounded in absolute value by a constant that converges 
to as e I 0. 

At this point, we would like to replace the currents W Xj i appearing in the first term of 
the right hand side of (3.15) by its convolution with uj. In order to do it we need some 
notation. For a local function g denote by g k the convolution 

{r x g) k (v) = j^-r E { r x+y9)(v) , 
1 fc| ?/eA fe 

when g(r)) = r)(0) — 1/2, we shall denote (r x g) k simply by i] k (x). Observe that, i] k (0) + 1/2 
is precisely the empirical density in A&. Furthermore, for e > 0, let A^ = {— e^ 1 + k, e^ 1 — 
k} x 71-°^ 1 and = A e \ A^. For x e X = Z x n^ 1 denote by W x j the current defined on 
X = {0, 1} X by 



W (n) - S Wx ^) if X,X + 6i £ A £ 

Wx,i{V)- | o if x orx + e* ^ A e 



with this notation we have 



E{ =e 



d-2 



I \ E E d £ ' l X(t,ex){r x W ,) k - £*$ 1 {f t - 9()(r,)<h, 

J [ i=l xgAJ J 



+£d ~ 2 1 EE^^W 

^ I »=i xerj 



(3.18) 

(/,-*f)fa)di/ 



:= E\ + E\ . 



Let us summarize what we have done so far: for any fixed A > and a > 0, we got the 
following bound for the entropy derivative 

j/~ 2s ^) < Ae d ~ 2 Sl (t) + E 2 + E\ + E\ 

- e d ~* {^o(v^f , *f ^1/2) + l -^D h {Jh u %u 1/2 ) } ( 3 - 19 ) 
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/From the fact that d £ ' l \{t 1 ex) has a compact support, Taylor expansion and (3.17), 
the term Ef can be written as 



+ / EE F' 1 ^ ™) (^?) fc - (ft - *n 

^ [ i=l xgAJ J 

+ £ (1). 



(3.20) 



Following the method of [EMY1] we now replace the currents (Wg^)* in the bulk by 

a linear combination of the gradients (i](ej) — i](0)) h , j = 1, • • • ,d. This requires some 

notation. For 16 A' and 1 < i < d, denote by the conditional expectation of W^J 
given the density of particles on A x 



M k x M = E 



where A Xjfc = {x + y : y £ A k }. One can compute M^^rj) easily. It is given by 
Furthermore, for e > 0, 1 < % < d and rj E X e , let 

i=i 

where Djj = Djj — 5jj and the diffusion matrix D i j is defined in Section 2. 

The next result is the main step towards the proof of the bounds of the entropy. 
Theorem 3.1. For all A > and any probability density f with respect to v, 



lim limlim<£^ / ( £ £ ^A(t, e^V^/fa)*/ - Ae^D^Jf ,v) 

n— >oo £— »oo £— >0 | / \ — * ^— ' / 

«=1 xGA* 



V < 0. 



This Theorem was proved in [EMY1], [LY]. The next lemma takes care of the terms 
close to the boundary (Ef) 

Lemma 3.2. For all A > and any probability density f with respect to u, 

hmli^ L d - 2 f ^' ix ^ £x K^*,i) k f d » ~ A£d ~ 4 My/f,v)\ <o. 
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The proof of this Lemma will be given in next section. 

Notice that, in the previous results, the Dirichlet form .D (\/7t> z/ ) appeared while, in the 
derivative of the entropy, we got D Q (^/h^ , tyf v). The next lemma allows us to replace the 

one by the other. Moreover we will also need some estimate about D (\J~$f , v). 
Lemma 3.3. There exist two constants C and C such that 



Proof. The second inequality follows by inspection. For two nearest neighbor sites x, y in 
A e , denote by S XtV the operator defined by 

{s s ,yf)(v) = f(v x ' y )-m- 

and write the Dirichlet form -Do(vTt) v) as 

d „ 

MVfw) = E E / ^^{(^^(^f)-')^)^^) 

i=l x,x+ei&K e 

+ ( y S x , x+ei ^hf)( V )} 2 ^( V )du. 
The elementary inequality (a + b) 2 < 2 (a 2 + 6 2 ) and a change of variables give that 
e^Doiy/fuv) <2e d - 4 D (yfhf^fu) 

d „ 

E E f<, 

i=l i.i+e.gAf 



□ 



Applying Lemma (3.2) with densities ft then \I/f and using Lemma 3.3, we obtain the 
following bound for Ef: for any A > 

£ x 4 < A^AKv^, + + r(e, k, A) (3.21) 

where 

lim \imr(e,k,A) < 0. 
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With the notations introduced before Theorem 3.1, a summation by parts and a Taylor 
expansion permit to rewrite the quantity Ef as 



+ E d - 2 



E\ = e^ 1 /{EE dl Uj \(t,ex)D i , j r 1 k (x)} (f t - *f) fa)di/ 

l<i,j<dxeAI 

7 { E E ^A(^)M^fa)} (/ t - tff) fa)*, 
^"7 { E E PtXfrexfeVM} (f t - *f) fa)*/ 
^7 {E E 9 B,i A(t,ex)(t,ex)i 7 *(x)ni(ex)}(/ t -*f) fa)di/ 



(3.22) 



+ R £ (t) + o £ (l) , 



where <9A^ stands for the boundary of A*, = (ni(ex), 0, • • • ,0) for the outward unit 

normal vector to e(dA*) at ex and i2 e (i) has been defined in (3.16). 

We first estimate the third line in the formula above. Using Theorem 3.1 and Lemma 
3.3, we get, for any A > 0, 



where 



£d ~ 2 j {EE ^X(t,ex)r x Yt(n)} (ft-*?) (r,)du< 
< Ae d - 4 D ( y/f t , u) + A(C' + C ) + r(e, k, n, A) 



lim lim \imr(e,k,n,A) < 0. 

n— >oo ^— >oo e^O 



(3.23) 



We now examine the second line in (3.22). We can write as 

5iQfc, e (a;) + 2e5iA(*, ex)r] k (x) + const + 0(k~ d ) + 0(£ 3 £; 2 ) , 

where 

Q k ,e(v, x ) = {v k ( x ) - e\{t,ex)) 2 

and const stands for a term independent of the configuration. By the entropy inequality, 
for any bounded function G:lx 7r d_1 — > R, for all q > 

£ d - 2 1 E c( £ x)g fe , £ fa^) (/* - *f) fa)*/ 

< g-V- 2 Sl (t) +g- 1 e d - 2 logf?*f [expg E ^fa^)] 
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with 

Qk,e(Vi x ) '■= Qk,e(V, X) - E** [Q k ,e(Vi x )] ■ 

^From large deviations estimate (Lemma 3.1. in [EMY1]) there exists go > 0, such that, for 
all q < g , 



lim lim£ d - 2 logE^[exp q V Q k , £ (r],x)} = 0. 

In conclusion, we got the following inequality for the second line in (3.22) 

£d ~ 2 /EE ^\{t,ex)MlM (/ t - ttf) 
i=i xeAJ 

^ £d_1 /EE ea; ) ( 2A (*» ^Kfr) (/* - *f ) 



»=1 zeAj 

+ g-V- 2 Si(t)+r £ifc (g), 

where, for all g < y, r £jk (q) converges to when e { and £ f oo. 

To deal with the fourth line in (3.22) (boundary term), we need the following lemma 

Lemma 3.4. Fix a bounded function G : R x 7r d_1 — > R and #i G [— + k, —e^ 1 + 2/c] U 
[e^ 1 — 2k, e" 1 — k]. Then there exists some constant C > such that for all A > 0, 



e d - 2 I J2 GMxiMrftx^y) (f t - ttf) ( V )du < 



The proof of this lemma is postponed to the next section. 

We now turn to the proof of Proposition 3.1 and we first consider the term E 2 . Since the 
function (d t X) is equal to outside Qe (cf. formula (3.5) and below), a simple computation 
shows that 

E 2 = -e^l d £Lf t dv 
J *f 

= ~ £d ~ 1 I J2 d ^ X *^ £x ^ (/t"*f) (V)du 

= - £d_1 / E ^ A (*> ^0*0 (/* - *f ) fa) ^ • 



x:sx£Qg 



To conclude the proof of proposition 3.1, we integrate the inequality (3.19) from to t. 
Combining with the above estimates, we obtain that there exist constants K\ > 0, K 2 > 0, 
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(3.24) 



HYDRODYNAMICS OF A DRIVEN LATTICE GAS 

K 3 > and c > such that, for all small enough A > 0, 

£ _2 si(£) < i^i£~ 2 / si(u)du 
Jo 

-K 2 e d - 4 J ds{D (y/}g,*?v)+D b (y/h a ,* 9 v)} 
+ e d - 1 fds{ V ff( a ,ex) / - tff )&/} 

+ AK 3 e d ~ 4 [ D ( y /f~s,i')ds + r(e,k,n,A) + cA, 
Jo 

where lim lim limr(e, k, n, A) < for all A > 0, and the function H is given by 

n^oo <!— »oc e~^0 

d 

H(s,ex) = ^25id ei (X(s,ex)) 2 + ^ Aj^ i) « j A(s,ex) - d s \(s,ex) . 

i=l l<i,j<d 

Recalling that \(t,u) satisfies (3.6) in fi 2 e and using the notation R £ defined in (3.16), we 
have 

e d ~ x J2 H ( s > £x ) j V k {x)Us-^> F s )dv = E d - 1 J2 H ^ £X ^ J V k (x)(f s -^)du+R £ (s). 



We will prove in the next section: 



Lemma 3.5. For any bounded function G defined on A and for any density f with respect 
to v, we have for any A > 



<Ae d ' 4 [D b {^h s ^ s u)+D Q {^T s ,u)] + Cd'\\G\\ ^ l|G||c 



+ 1 



for some positive constant C . 

Therefore, noticing that ||#||oo < const 9~ 2 , we get for any small enough A > 



£ 2 si{t) < K\S 2 / si(u)du 
Jo 



K 2 e a -* / ds 



J^ds{D ( y/hf, ttf u) + D b ( y/h„ } 

+ AX 3 £ d - 4 jf rfs A) ( v 7 ^ ^) + r(e, A;, n, A) + c (a + 
where the constants previously defined may have changed their values. 



(3.25) 
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To conclude, we use Lemma 3.3 to control the Dirichlet form e d 4: D (^/J s ~,h') with 
E d ~ 4 Do( y/hfj', tyfv)- The error obtained from this replacement is a constant but we notice 
that there is the small factor A in front of this term. Therefore there exists a constant 
d > such that for any small enough A > 0, 

£ _2 si(£) < KiE~ 2 / si(u)du 
Jo 

-K 2 e d ~ 4 J ds {D (y/h£,*?u)+D b (y/h a ,y a u)} 

(q 
A+ A 

then we choose A = A(9) vanishing with 9 in such a way that \im9A~ 1 (9) = (e.g. 

A{9) = V9). Finally we consider the successive limits e — » 0, £ oo, n — > oo, 9^0 and 
we apply Gronwall lemma. 

□ 




4. Estimates on boundary terms 



Proof of Lemma 3.2. From the definition of W X j we have that the term Ef in (3.18) 
can be written as 

^ / 1 E E H£,i ^ £x ^jh E \ (/* - *0(^ 

l<=ixer* 1 fc| s/eAfcW J 

where A fe (x) is the block of rectangular shape which is the set of y G A k (x) (cube centered 
in x) such that y and y + e\ belong to A e . Let M xi (rj) be the conditional expectation 
of W x ^,x G A £ given the density of particles on Afc(x). Repeating the argument given in 
the paragraph following inequality (3.23), it suffices to prove that, for all A > 0, for any 
bounded function J and any probability density / with respect to u, 



lim lim < e 

n^oo e— >0 



d-2 



/d ^ ^ 

( E E J ^ £x ) [i W ^ * ^) W ~ ) fdy 

-As d - 4 D (yff)\=0 



We need the following definition: 

Let fj,£ be the canonical measure in the block A^ with given density rf. For any {x, y} C 
A^, we introduce the Dirichlet form Dq{K) as 



D^ y \h) = j r {xM (rj)(h(rj^) - h( V )) 2 dpt e 
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and we define the finite volume variance 

V t {Grf) = (21 + l)- d (J2(r x G - EfflrflX-CP)- 1 £(r x G - E[G\ V %, r 

\x\<e \x\<e 

where is the symmetric part of the generator £ £j0 restricted to the box A e . The proof 
of Theorem 4.6 in [EMY1] is based on the following result: 

Lemma 4.1. For any cylinder function h there exist a constant C and a function C(h,£) 
vanishing for £ — > oo such that for any positive A and I = £ d+2 



I 



uj * [r x h - E[h\ V k (x)]f t dLi k -A—- V Di v ' z} (y/f) 

{2/,2>CAfc(a;) 



< V l {h,rf)f t dLik + E 2 C{h,e) 



The proof of this Lemma is given in [EMY1] for square blocks but it can be extended 
easily to a rectangular shape provided that the volume of the block is of order k d . Using 
Lemma 4.1 we prove Lemma 3.2 by taking the expectation with respect to u, multiplying 
by e d ~ 2 J, summing over x G and noting that the number of terms in T k is e~ d+1 k. 

Proof of Lemma 3.4. Fix e" 1 — 2k < x\ < e^ 1 — k and write the sum as 

e d ~ 2 f G{e(x 1 ,y))i] k (x 1 ,y)(f t - V()dv 

J r- d-1 

= £ T"E / E {G(e(x u y))) k r,(x 1 + z 1 ,y)(f t -*[)*,, 

\zi\<k J yen d-l 

where (G(e(xi,y))} k is the (d — l)-dimensional convolution in the variable y E 7Tg -1 . 
Therefore, to prove the lemma, it is enough to show that, if H : R — > R is a bounded 
function, then for every e^ 1 — 3k < z\ < e^ 1 , and A > 

e d - 2 ( H(e(z 1 ,y))r j (z 1 ,y)(f t -^)du 

< e d - 4 ^ (D b {^h u %u) + D (y/f t , !/)) + const eh]^ + lj , 

Let e' 1 — 3k<z 1 < e' 1 and decompose the left hand side of the last inequality into two 
terms Bi and B 2 

d ' 2 S E H ^ v)) y) - it 6 ' 1 * v)) (/* - 

J ,_ d-1 

I H(e( Zl ,y)) V (e-\y)Ut-^t)dv 

, r- d-1 



s- CI 

y€ir e 
d-2 



B ± + B 2 
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The term B 2 is the simplest one. From an integration by parts and Schwartz inequality, 
it is bounded, for all A > 0, by 

^<^ 2 E {^-^-^(^M+e^llffllSoj, (4.26) 

d-l V ' 

where for a function /, a positive measure \i and x G T e 

D b x (f,n) = ((f(v x )-m) 2 ),. 

We now consider B\. By Schwartz inequality 

fil = H H « z ^ y)) + * + l,v)-v(t + z, y)) (ft - *f )*/ 

< e d ~ 2 |i^D (v^, + -^rW H Wl\ + ekC \\H\U 

for some constant C > 0, where we have used in the last inequality, for the second term 
corresponding to the integration with respect to \l/f u, integration by parts and Taylor 
expansion. □ 

Proof of Lemma 3.5. The summation over the set {x , ex G Q £ k \ can be divided in 
two similar terms. We consider the one where the first coordinate is such that e~ 1 (l — 29) < 
Xi < e^ 1 — k, the second term is handled in the same way. If we repeat the arguments 
used in the proof of the Lemma 3.4, we obtain for £ _1 (1 — 29) < x x < e' 1 — k, 



e d - 



1 G « x ^y)) j v\xi,y){m-^t)dv 



yen, 

<zU d -^ e (D (y/f u u) + D b (Vh t ,%u)) +C7^G|U(^ + l)) 
To conclude the proof of the lemma, we just have to take the sum over x\. 



5. Comments 

We conclude with a few generalizations: 
1) Assumption 8± > 0. 

If 5i — we need to introduce a different boundary generator to fix the density on the 
boundary. We make the choice which is usually done for the symmetric case, a death and 
birth process Cb acting on each site of the boundary such that it is reversible with respect 
to the one site measure with density p(ex) = 1/2 + eb(ex) 

Ufiji) = £ [ P (ex)(l - r) x ) [/(rf '+) - f{ v j\ + (1 - p(ex))r) x [/(^") - f(r,)} ' 

where i] x,+ and i] x, ~ are the configurations obtained from 77 by adding or removing a particle 
at site x. 



HYDRODYNAMICS OF A DRIVEN LATTICE GAS 



21 



2) General domain. 

We generalize now the model to the case of a macroscopic system in a smooth bounded 
convex domain of M. d . We introduce the boundary generator acting on the boundary T £ 
as in Section 2 with the boundary rates C b (ex,r]) chosen in such a way that a particle is 
added at the site x on the boundary when the site is empty, with rate \S ■ n\(l/2 + eb(ex)) 
if S ■ n, the scalar product of the vector 5 and the outward normal in x, is positive and 
removed when the site is occupied, with rate \S ■ n\(l/2 — eb(ex)) if the scalar product is 
negative. For x G T £ 

C\ex, rj) = \5 ■ n\ + e b(ex))(l - rj(x)) 1{6 ■ n(ex) > 0} 

+ \5 ■ n\ (- - e b(ex))i](x) 1{5 ■ n(ex) < 0}. 

These rates fix the value of the particle density on the macroscopic boundary to be 1/2 + 
e b(x). This choice is sufficient to remove the entropy flow generated by the drift of ASEP. 
Since S-n(ex) can be zero in some points for a general domain we have to add, on the basis 
of the previous remark, also a reversible boundary generator to fix the density in these 
sites. 

3) Different boundary conditions. 

To remove the entropy flux it is possible to choose a generator C b which is reversible instead 
of using the non reversible generator C Et b, but then we have to speed up it by a factor s~ 3 
with respect to the jump process. The total generator speeded up is then 

e 2 C £ = e 2 C £: o + ^ 5 £b- 

We note that the bad boundary terms in our proof are eliminated by means of a cancella- 
tion. If we use only the reversible generator e~ 2 Cb on the boundary they will be controlled 
by the Dirichlet form associated to £&, by a generalization of the Lemma 3.3. 
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